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SPIN-DOWN MEASUREMENT OF PSR J1852+0040 IN KESTEVEN 79: 
CENTRAL COMPACT OBJECTS AS ANTI-MAGNETARS 
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ABSTRACT 

Using XMM-Newton and Chandra, we achieved phase-connected timing of the 105 ms X-ray pulsar 
PSR J1852+0040 that provides the first measurement of the spin-down rate of a member of the class 
of Central Compact Objects (CCOs) in supernova remnants. We measure P = (8.68 ± 0.09) x 10 -18 , 
and find no evidence for timing noise or variations in X-ray flux over 4.8 yr. In the dipole spin-down 
formalism, this implies a surface magnetic field strength B s = 3.1 x 10 10 G, the smallest ever measured 
for a young neutron star, and consistent with being a fossil field. In combination with upper limits on 
B s from other CCO pulsars, this is strong evidence in favor of the "anti-magnetar" explanation for 
their low luminosity and lack of magnetospheric activity or synchrotron nebulae. While this dipole field 
is small, it can prevent accretion of sufficient fall-back material so that the observed X-ray luminosity 
of L x = 5.3 x 10 33 (<i/7.1 kpc) 2 erg s _1 must instead be residual cooling. The spin-down luminosity 
of PSR J1852+0040, E = 3.0 x 10 32 erg s~\ is an order-of-magnitude smaller than L x . Fitting of 
the X-ray spectrum to two blackbodies finds small emitting radii, R\ — 1.9 km and R2 — 0.45 km, 
for components of kT\ = 0.30 keV and kT^ = 0.52 keV, respectively. Such small, hot regions are 
ubiquitous among CCOs, and are not yet understood in the context of the anti-magnetar picture 
because anisotropic surface temperature is usually attributed to the effects of strong magnetic fields. 
Subject headings: ISM: individual (Kes 79) — pulsars: individual (PSR J0821-4300, IE 1207.4-5209, 
PSR J1852+0040) — stars: neutron 



1. INTRODUCTION 

The class of relatively faint X-ray sources known as 
central compact objects (CCOs) in supernova remnants 
(SNRs) are apparently isolated neutron stars, which we 
define here by their steady flux, predominantly ther- 
mal X-ray emission, lack of optical or radio counter- 
parts, and abse nce of a surroundin g pu lsar wind nebul a 
(see reviews by iPavlov et all 120041 and iDe Luca1l2008f) . 
Table [T] lists basic data on the well-studied CCOs, as 
well as candidates whose qualifications are not well es- 
tablished. Of the seven most secure members, three 
are defi nitely pulsars with p e riods of 0.105, 0-112, and 
0.424 s (IGotthelf et alj|200a IGotthelf fc Halper"3 l2009t 
iZavlin et all I2000D . Until now, no spin-down was de- 
tected from a CCO, which is most simply interpreted 
as indicating a weak surface dipole magnetic field B s . 
In two cases, PSR J1852+0040 (P = 0.105 s) and 
IE 1207.4-5209 (P = 0.424 s), the upper limits on 
B s from spin period measuremen ts are 1.5 x 10 11 G 
and 3.3 x 10 11 G resp ectively (jHalpern et all 120071 : 
IGotthelf fe Ha lpcrn 2007), smaller than that of any other 
young neutron star. We also found B s < 9.8 x 10 11 G 
from a pair of observatio ns of PSR J0821— 4300, t he 
0.112 s pulsar in Puppis A (jGotthelf fc Halpern|[2009l ). 

Important implications of these results are that the 
birth periods of the CCO pulsars are not significantly 
different from their present values, and that their spin- 
down luminosities are, and have always been, insufficient 
to generate significant non-thermal magnetospheric emis- 
sion or synchrotron nebulae. A corollary is that the so- 
called "characteristic age" r c = P/2P that is used to 
approximate the true age of pulsar has, for CCOs, no 
meaning, being at least millions of years for pulsars that 
are demonstrably in supernova remnants that are only 



a few thousand years old. It is also reasonable to sup- 
pose that some isolated radio pulsars with weak magnetic 
fields that have characteristic ages of millions of years 
may actually be former CCOs that are only moderately 
aged. 

Most of the properties of the CCOs c an thus be ex- 
plained by an "anti-magnetar " model (jHalpern et al.l 
120071: IGotthelf fc Halperif[2008l ). including the possibil- 
ity that their weak magnetic fields are causally related to 
their slo w rotation periods at birth th rough the turbulent 
dynamo (|Thompson fc Duncanlll993l ) that generates the 
magnetic field. ( See iSpruitl 12008 for a review of possible 
mechanisms for the origin of magnetic fields in neutron 
stars.) While CCOs are inconspicuous relative to ordi- 
nary young pulsars and magnetars, the fact that they are 
found in SNRs in comparable numbers to other classes 
of neutron stars implies that they must represent a sig- 
nificant fraction of neutron star births. Considering that 
the CCO in Cassiopeia A is the youngest known neutron 
star (330 yr), and postulating that a CCO in SN 1987A 
explains why its pulsar has not yet been detected, argues 
that anti-magnetars are potentially a populous class. 

The compact X-ray source CXOU J185238.6+004020 
was discovered in the center of the SNR Kes 79 by 
ISeward et alj (|2003D . In previous papers, we reported 
the discovery of 105 ms p ulsations from that CCO, now 
named PSR J1852+0040 (|Gotthelf et al.ll2005l . Paper 1), 
and the first few observations that established only an 
upper limit on its period derivative c orresponding to 
B s < 1.5 x 10 11 G (Ha lpern et ail 120071 Paper 2). Here, 
we present a dedicated series of timing observations of 
PSR J1852+0040 that constitute the first definite mea- 
surement of the spin-down rate of a CCO pulsar, confirm- 
ing its unusually small dipole magnetic field and support- 
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TABLE 1 

Central Compact Objects in Supernova Remnants 



ceo 


SNR 


Age 


d 


P 




B s 




References 






(kyr) 


(kpc) 


(s) 


(%) 


(10 11 G) 


(erg s L ) 




RX J0822.0-4300 


Puppis A 


3.7 


2.2 


0.112 


11 


< 9.8 


6.5 x 10 33 


1,2 


CXOU J085201. 4-461753 


G266. 1-1.2 


1 


1 




< 7 




2.5 x 10 32 


3,4,5,6,7 


IE 1207.4-5209 


PKS 1209-51/52 


7 


2.2 


0.424 


9 


< 3.3 


2.5 x 10 33 


8,9,10,11,12 


CXOU J160103. 1-513353 


G330.2+1.0 


> 3 


5 




< 40 




1.5 x 10 33 


13,14 


1WGA J1713.4-3949 


G347.3-0.5 


1.6 


1.3 




< 7 




~ 1 x 10 33 


7,15,16 


CXOU J185238. 6+004020 


Kes 79 


7 


7 


0.105 


64 


0.31 


5.3 x 10 33 


17,18,19,20 


CXOU J232327.9+584842 


Cas A 


0.33 


3.4 




< 12 




4.7 x 10 33 


20,21,22,23,24 


XMMU J172054.5-372652 


G350.1-0.3 


0.9 


4.5 








3.4 x 10 33 


25 


XMMU J173203.3-344518 


G353.6-0.7 


~ 27 


3.2 








1.0 x 10 34 


26,27,28 


CXOU J181852.0-150213 


G15.9+0.2 


1-3 


(8.5) 








~ 1 x 10 33 


29 



References. — (1 ) IHui fc Becker M2006t (2) IGotthelf fc Halpernll200a (3) ISlane et al.N2001l; (4) [Kargaltsev et al.M2002t (5) 
Bamba et al.l2005l: (6lllvudin et al.l200a miDe Luca 120081: fS l lZavlin et al.120001: (9 1 IMereghetti et al.l2002l : (10 1 | Bjgnami et all 
2003; (11) De Luca ct al. 2004; (12) Gotthelf & Halpern 2007; (13) Park et al. 2006; (14) Park et al. 2009; (15) Lazcndic et al. 
2003; (16) Cassam-Chcna'f ct al. 2004; (17) Seward ct al. 2003: (18) Gotthelf ct al. 2005; (19) Halpern ct al. 2007; (20) this paper: 
(21) Pavlov et al.H200Ct (22 ) I ChaJtrabaxtv et alJI200H ( 23) IMereghetti et al. l po'oa (24) [Pavlov fc Luna| | 2009| ; (25) IGaensler et all 
12008 ; (26) ITian et alJI2008t (27) IAcero et al.ll2009t (28) lHalpern fc Gotthelfcooa (29) lRevnolds et aljhood . 

Note. — Above the line are seven well-established CCOs. Below the line arc three candidates. 
a Upper limits on pulsed fraction are for a search down to P — 12 ms or smaller. 



ing the anti-magnetar model. The plan that was devised 
to obtain the needed series of time-constrained observa- 
tions is described in Section 2. Analysis and results of 
the timing data are presented in Section 2.1, and the 
high-quality X-ray spectrum and pulse profile that were 
obtained from the summed observations are the subject 
of Section 2.2. Discussion of the anti-magnetar model, 
and possible explanations for the X-ray spectrum appear 
in Section 3, followed by the conclusions in Section 4. 

In this paper, we adopt a distance of 7. 1 kpc to Kes 79 
from H I and OH absorptio n studies (|Frail fc Cliftonl 
119891 : [Green fc Dewdnev 1992) u pdated using the G alac- 
tic rotation curve of ICase fc Bhattacharval (|T998V We 
also assume the dynami cal estimate of 5 .4-7.5 kyr for 
the age of the SNR from lSun et all (|200l . 

2. X-RAY OBSERVATIONS 

After the first four timing observations of PSR 
J1852+0040 in 2004 and 2006 (Paper 2), it was evident 
that the frequency derivative / was too small to mea- 
sure without obtaining a phase-coherent series of obser- 
vations spanning several years. The most economical way 
to measure / is to begin with a logarithmically spaced se- 
quence of observations that maintains the absolute cycle 
count <j){t) /2n and measures the frequency / with in- 
creasing accuracy, until the small quadratic term in the 
phase ephemeris 

f 1 = f(t-t )+ 1 -f(t-t ) 2 + ... 
Zir I 

begins to make a significant contribution to the phase. 
In this case, it was also deemed possible (Paper 2) that 
accretion from fall-back material could contribute tim- 
ing irregularities known as torque noise, of a magnitude 
comparable to the existing upper limits on dipole spin- 
down. This made it all the more important to obtain a 
well-sampled ephemeris that could test for such effects. 

Unfortunately, maintaining cycle count would require 
more observations classified as time-constrained than the 
Chandra X-ray Observatory allocates to any one project, 



while for XMM-Newton, semiannual visibility windows 
for this source are only 40 days long, not wide enough 
to securely bridge over the intervening 5 month gaps. 
It took until 2008 to implement a strategy that uses 
the two satellites in a coordinated fashion, filling two 
XMM-Newton visibility windows with six observations 
each of variable spacing, while requesting pairs of Chan- 
dra observations to bridge the gaps between and af- 
ter the two XMM-Newton windows. In this manner, a 
phase-coherent timing solution spanning 1.7 yr could be 
achieved, which could also be extrapolated backward to 
incorporate the earlier observations. All but one of the 
approved observations in 2008-2009 were performed as 
planned. Due to a loss of contact with the XMM-Newton 
spacecraft, the last observation of 2008 was rescheduled 
to 2009, but this did not compromise the success of the 
program. We were thus able to obtain a fully coher- 
ent timing solution incorporating all of the observations 
listed in Table [2j including the earliest ones, spanning 
4.8 yr in total. 

All of the XMM-Newton observations used the pn de- 
tector of the European Photon Imaging Camera (EPIC- 
pn) in "small window" (SW) mode to achieve 5.7 ms 
time resolution, and an absolute uncertainty of « 3 ms 
on the arrival time of any photon. We processed all EPIC 
data using the emchain and epchain scripts under Science 
Analysis System (SAS) version xmmsas_20060628_1801- 
7.0.0. The leap second at the end of 2008 was inserted 
manually. Simultaneous data were acquired with the 
EPIC MOS cameras, operated in "full frame" mode. Al- 
though not useful for timing purposes because of the 2.7 s 
readout, the location of the source at the center of the 
on-axis MOS CCDs allows a better background measure- 
ment to test for flux variability, an important indicator 
of accretion, than the EPIC-pn SW mode. 

The Chandra observations used the Advanced Cam- 
era for Imaging and Spectroscopy (ACIS) in continuous- 
clocking (CC) mode to provide time resolution of 
2.85 ms. This study uses data processed by the pipeline 
software revisions v7.6.9-v8.0. Reduction and analysis 
used the standard software package CIAO (v3.4) and 
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TABLE 2 

Log of X-ray Timing Observations of PSR J1852+0040 



Mission 


Instr/Mode 


ObsID/Seq# 


Date 




Exposure 


Start Epoch 




Period a 


O 










CUT) 




(to) 


(MJD) 




(ms) 




XMM 


EPIC-pn/SW 


0204970201 


2004 


Oct 


18 


30 


g 


53296.001 


104 


912638(39) 


121 


7 


XMM 


EPIC-pn/SW 


0204970301 


2004 


Oct 


23 


30 


5 


53301.984 


104 


912612(55) 


77 


2 


XMM 


EPIC-pn/SW 


0400390201 


2006 


Oct 


08 


29 


7 


54016.245 


104 


912610(47) 


92 


4 


Chandra 


ACIS-S3/CC 


6676/500630 


2006 


Nov 


23 


32 


2 


54062.256 


104 


912592(40) 


94 





XMM 


EPIC-pn/SW 


0400390301 


2007 


Mar 


20 


30 


5 


54179.878 


104 


912600(43) 


107 


5 


Chandra 


ACIS-S3/CC 


9101/500964 


2007 


Nov 


12 


33 


I 


54426.674 


104 


912615(40) 


95 


8 


Chandra 


ACIS-S3/CC 


9102/500965 






16 


■ >i 


2 


54628.080 


104 


912593(40) 


115 


5 


XMM 


EPIC-pn/SW 


0550670201 


2008 


Sep 


19 


21 


2 


54728.758 


104 


912563(84) 


61 


7 


XMM 


EPIC-pn/SW 


0550670301 


2008 


Sep 


21 


31 





54730.066 


104 


912594(46) 


84 


2 


XMM 


EPIC-pn/SW 


0550670401 


2008 


Sep 


23 


34 


8 


54732.079 


104 


912573(42) 


80 


3 


XMM 


EPIC-pn/SW 


0550670501 


2008 


Sep 


29 


33 





54738.016 


104 


912596(40) 


97 


4 


\f 71 IT 71 IT 

XMM 


EPIC-pn/SW 


0550670601 


2008 


Oct 


10 


36 





54750.006 


104 


912641(40) 


90 





Chandra 


ACIS-S3/CC 


9823/501015 


2008 


Nov 


21 


30 


1 


54791.803 


104 


912645(55) 


69 


8 


Chandra 


ACIS-S3/CC 


9824/501016 


2009 


Feb 


20 


29 


6 


54882.483 


104 


912591(53) 


80 


1 


XMM 


EPIC-pn/SW 


0550671001 


2009 


Mar 


16 


27 





54906.255 


104 


912610(48) 


97 


9 


XMM 


EPIC-pn/SW 


0550670901 


2009 


Mar 


17 


26 





54907.609 


104 


912592(60) 


82 


9 


XMM 


EPIC-pn/SW 


0550671201 


2009 


Mar 


23 


27 


3 


54913.582 


104 


912616(46) 


104 


8 


XMM 


EPIC-pn/SW 


0550671101 


2009 


Mar 


25 


19 


9 


54915.654 


104 


912544(85) 


75 


8 


XMM 


EPIC-pn/SW 


0550671301 


2009 


Apr 


04 


26 





54925.539 


104 


912590(52) 


97 


9 


XMM 


EPIC-pn/SW 


0550671901 


2009 


Apr 


10 


30 


5 


54931.535 


104 


912620(48) 


82 


9 


XMM 


EPIC-pn/SW 


0550671801 


2009 


Apr 


22 


28 





54943.897 


104 


91262l(69) 


67 


i 


Chandra 


ACIS-S3/CC 


10128/501055 


2009 


Jun 


02 


33 


2 


54984.888 


104 


912596(37) 


117 


9 


Chandra 


ACIS-S3/CC 


10129/501056 


2009 


Jul 


29 


32 


2 


55041.231 


104 


912633(40) 


105 


6 



a Barycentric period derived from a test. The Leahy ci al. (1983|) uncertainty on the last digits is in parentheses. 



CALDB (v3.4.2). The photon arrival times in CC mode 
are adjusted in the standard processing to account for 
the known position of the pulsar, spacecraft dither, and 
SIM offset. Absolute accuracy of the time assignment in 
Chandra CC-mode is limited by the uncertainty in the 
position of the pulsar, which was determined in Paper 2 
from an ACIS image. The typical accuracy of w 1 pixel 
then corresponds to an uncertainty of w 3 ms, which is 
similar to the XMM-Newton accuracy, and is 0.03 rota- 
tions in the case of PSR J1852+0040. We will show that 
the measured dispersion in pulse arrival times is compa- 
rable to this uncertainty. 

2.1. Results of Timing Analysis 

For each observation in Table we transformed the 
photon arrival times to Barycentric Dynamical Time 
(TDB) using the pulsar coordinates determined in Pa- 
per 2 and reproduced in Table [3] Diffuse SNR emission 
is a significant source of background. To maximize the 
pulsar signal-to-noise ratio, we used a source extraction 
radius of 12" for the XMM-Newton observations, and five 
columns (2."4) for the Chandra CC-mode data. An en- 
ergy cut of 1—5 keV was found to maximize pulsed power. 
The pulse profile and value of the period i n each observa- 
tion w as derived from a Z\ periodogram (jBuccheri et alJ 
1983), a choice of harmonics that was found to mini- 
mize the uncertainties. The resulting profiles were cross- 
correlated, shifted, and summed to create a master pulse 
profile template. Individual profiles were then cross cor- 
related with the template to determine the time of arrival 
(TOA) at each epoch. 

Starting with the dense series of XMM-Newton ob- 
servations in 2008 September, the TOAs were itera- 
tively fitted using the TEMPO 1 software. We fitted 
the three observations from September 19-23 to a lin- 

1 http:/ /www. atnf.csiro.au/research/pulsar/tcmpo 



ear ephemeris, and added TOAs one at a time using a 
quadratic ephemeris, finding that the new TOA would 
match to < 0.1 cycles the predicted phase derived from 
the previous set. After adding the final observation to the 
ephemeris, we then worked backward in time from 2008 
September to 2004 October until all 23 observations had 
been included. The quadratic term contributes —8.7 cy- 
cles of rotation over the 4.8 yr span of the ephemeris, 
which yields the small uncertainty in P discussed be- 
low. The phases and errors were determined by cross- 
correlating with a final iterated template of co-added 
profiles produced from the best fitting ephemeris given in 
Table [3] Figure [1] shows the residuals of the individual 
observations from the best fit, which demonstrates the 
validity of the solution. The weighted rms of the phase 
residuals is 3.4 ms, or 0.032 pulse cycles, which is compa- 



PSR J1 852+0040 TIMING SOLUTION RESIDUALS 
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Fig. 1. — Pulse phase residuals for the 23 timing observations of 
PSR J1852+0040 from the best fitting ephemeris of Table|3] Error 
bars arc 1 a. 
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TABLE 3 

Spin Parameters of PSR J1852+0040 



Parameter 



Value 



Right ascension, R.A. (J2000) a 
Declination, Decl. (J2000) a 
Epoch (MJD TDB) b 
Spin period, P 
Period derivative, P 
Valid range of dates (MJD) 
Surface dipole magnetic field, B a 
Spin-down luminosity, E 
Characteristic age, r c 



18 h 52 m 38 s .57 

+00°40'19".8 

54597.00000046 

0.104912611147(4) s 

(8.68 ±0.09) x 10- 18 

53296 - 55041 

3.1 x 10 10 G 

3.0 X 10 32 erg s" 1 

192 Myr 



a Measured from Chandra ACIS-I ObsID 1982 (Paper 1). Typ- 
ical Chandra ACIS coordinate uncertainty is 0. 6. 

Epoch of ephemeris corresponds to phase zero in Figure [2] 

rable to the individual measurement errors. There is no 
evidence of timing noise or higher derivatives in the resid- 
uals. The variation of Zf(pulsed power) among the ob- 
servations listed in Table[2]is also consistent with statisti- 
cal expectations. Figure [5] shows the summed pulse pro- 
file from the f 6 XMM-Newton observations, for which it 
is possible to measure a reasonably accurate background 
(unlike the Chandra CC-mode data). 

The sensitivity of the results to P and B s can be 
estimated analytically. If a pulsar is spinning down 
smoothly, then a coherent timing solution spanning time 
T will have an uncertainty in frequency of Sf = SP/ P 2 ~ 
0.1/T, and will be sensitive to a period derivative of 
P min = SP w 0.2(P/T) 2 . This in turn will measure 

B s = 3.2 x 10 19 V~PP~ with fractional uncertainty 



= 0.05 



SB, 
B s 



10 10 G 



i.o x io 3 



T 



pa 



B 2 T 2 



4.8 yr 



P 



0.105 s 



A coherent ephemeris for PSR J1852+0040 spanning 
4.8 yr has an uncertainty on P of SP « 1 x 10 -19 , 
corresponding to a 5cr detection limit of B s = 5 x 
10 9 G. In comparison, the measured values P = 
0.104912611147(4) s and P = 8.68(9) x 10" 18 have fitted 
uncertainties consistent with these analytic estimates, 
and imply in the dipole spindown formalism a surface 
magnetic field strength B s =3.1 x 10 10 G, a spin-down 
luminosity E = -IVlVL = Att 2 IP/P 3 = 3.0 x 10 32 erg s -1 , 

and characteristic age r c = P/2P = 192 Myr, all with 
negligible statistical error. 

2.2. Flux and Spectral Analysis 

In Papers 1 and 2 we presented several observations of 
PSR J1852+0040 that were consistent with steady flux 
and a blackbody spectrum of kT m 0.46 keV. The lumi- 
nosity of 3.4 x 10 33 d 2 1 erg s _1 corresponded to a black- 
body radius of only f» 0.8^7.1 km. This is typical result 
for a CCO, and indicates a concentrated hot spot that 
remains to be understood. Using the methods described 
in Paper 1, we now combine all 16 XMM-Newton obser- 
vations of PSR J1852+0040 into one spectrum in order to 
search for deviations from a blackbody that might indi- 
cate temperature variations on the surface, or cyclotron 



PSR J1 852+0040 



1.00 < E(keV) < 5.00 




TOTAL BACKGROUND 



0.5 1 1.5 2 

PHASE (MJD 54597.00000046) 

Fig. 2. — Summed pulse profile from the 16 XMM-Newton ob- 
servations of PSR J1852+0040 listed in Table |2] folded using the 
ephemeris of Table \3\ Phase zero corresponds to the TDB epoch 
in Table [3] The pulsed fraction after correcting for background is 
f p = 64 ± 2%. 

lines. EPIC pn and MOS spectra were fitted jointly. The 
accumulated exposure allows us to fit the spectrum over 
0.7 — 7 keV, an increase in coverage over the 1 — 5 keV 
range used in the previous papers. While a single black- 
body can still be fitted with a temperature of 0.46 keV, 
the high signal-to-noise ratio and expanded energy range 
of the summed spectrum reveal systematic deviations 
from the fit; the reduced xt — 1-38 for 240 degrees of free- 
dom is unacceptable (see Figure [3] and Tabled]). Similar 
to the results from other CCOs, we find that a fit to two 
blackbodies is significantly improved (xl — 1-07), with 
temperatures of kTi = 0.30 keV and kT 2 = 0.52 keV, 
while the corresponding radii R\ = 1.9^7.1 km and 
i?2 = 0.45 dj.i km still cover only a small fraction of the 
surface. Here we define radius using the phase-averaged 
luminosity Li :2 = 47ri? 2 2 crT^ 2 regardless of the unknown 
emission geometry, which could be, for example, a con- 
centric annulus. Adoption of the two-blackbody fit re- 
sults in a higher bolometric luminosity than the single 
blackbody, 5.3 x 10 33 d 2 A erg s" 1 vs. 3.0 x 10 33 d 2 7 A 
erg s _1 . We don't combine the Chandra spectra here for 
an independent fit because of the increased background 
and other uncertainties involved in analyzing CC-mode 
data for this faint source. 

Additional spectral models that were explored include 
the nonmagnetic n eutron star hydrogen atmosphere of 
iZavlin et al.l (|1996l ) , a nd a simplified Compt onized black- 
body, as described bv lHalpern et al.l (|2008f) for applica- 
tion to anomalous X-ray pulsars. Each of these models 
fits nearly as well as the two-blackbody model because 
they can eliminate the residual excess at high energy that 
is left by a single blackbody fit (see Figure [4] and Ta- 
bled]). For the neutron star atmosphere (NSA), we fixed 
the parameters M — 1.4 Mq and R°° = 13.06 km (cor- 
responding to R = 10 km), treating the effective tem- 
perature and distance as free parameters. As expected, 
the fitted temperature is smaller than those of the black- 
body models, but the fitted distance of 23.7 kpc is a 
factor of 3.3 higher than the known value. This indi- 
cates that in the atmosphere model only ~ 10% of the 
surface is emitting. We are not motivated to search for 
an atmosphere model that fits to the full surface area, 
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TABLE 4 

XMM-Newton Spectral Fits to PSR J1852+0040 



Parameter 


1313 


BB+BB 


NSA 


Comp 1313 


iV H (10 22 cm" 2 ) 


1.32 ±0.05 


l-»2_ lg 


L - 01 -0.07 


, „+0.09 


kTi (keV) 


0.460 ± 0.007 


0.30 ± 0.05 


0.287 ±0.007 


0.396 ± 0.014 


fll (km) 


0.72 d r .i 


1.9 d 7 i 


13.06 a 


0.93 d 7 .i 


kT 2 (keV) 




0.52 ± 0.03 






R 2 (km) 




0.45 dr.i 






r 






„ 7+4.0 
zo. / _ 2 g 


4.80+°;^ 


d (kpc) 








F a (0.5- 10 keV) b 


2.0 x 10~ 13 


2.0 x 10~ 13 


2.0 x 10~ 13 


2.0 x 10~ 13 


Lbol (erg s" 1 ) 


3.0 x 10 33 df. j 


5.3 x 10 33 d 2 x 




3.6 x 10 33 4 i 


x 2 H 


1.38(204) ' 


1.07(201) ' 


1.10(204) 


1.12(203) ' 



NOTE. — Spectrum fitted over 0.7 — 7 keV. Errors are 90% confidence for one interesting 
parameter (\ 2 — Xmin 2-7). 

a Fixed parameter — 13.06 km corresponding to R = 10 km with M — 1.4 Mq. 

" Absorbed flux in units of erg cm -2 s _1 . Average of EPIC pn and MOS. 

c Unabsorbcd, bolometric luminosity assuming d — 7.1 kpc. Average of EPIC pn and MOS. 



as was done for the Cas A CCO bv lHo fc Heinkel (|2009f ). 
because the larged pulsed fraction of PSR J1852±0040 
clearly requires a small hot spot. The same interpreta- 
tion attaches to the Comptonizcd blackbody model, in 



PSR J1 852+0040 
16 XMM-NEWTON OBS. 
327/466 ks EPN/MOS _ 





ENERGY (KEV) 



PSR J1852+0040 
16 XMM-NEWTON OBS. 
327/466 ks EPN/MOS 




ENERGY (KEV) 

Fig. 3. — Spectrum of the 16 summed XMM-Newton observa- 
tions of PSR J1852±0040. The EPIC pn spectrum is black, while 
the average of the MOS1 and MOS2 spectra is red. Top: Fit to 
a single blackbody model. Bottom: Fit to a double blackbody 
model. The parameters of the fits are given in Table [4] 



which the inferred blackbody radius is only 0.93 dj,i km. 

The large intrinsic pulsed fraction, f p = 64 ± 2% 
(defined as the fraction of counts above the minimum 
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Fig. 4. — Spectrum of the 16 summed XMM-Newton observa- 
tions of PSR J1852±0040. The EPIC pn spectrum is black, while 
the average of the MOS1 and MOS2 spectra is red. Top: Fit to 
a nonmagnetic neutron star hydrogen atmosphere model. Bottom: 
Fit to a Comptonized blackbody model. The parameters of the fits 
are given in Table [4] 
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Fig. 5. — Energy-resolved pulse profiles from the summed 16 
XMM-Newton observations of PSR J1852+0040 of Figure[2]decom- 
posed into four energy bands containing equal numbers of photons. 
Background has been subtracted. The pulse shape does not differ 
significantly as a function of energy. 

in the light-curve, corrected for background; see Fig- 
ure [2]) confirms that the emitting area must be small 
and far from the rotation axis, while the line of sight 
also makes a large angle to the rotation axis. The 
pulse shape does not appear to be a function of en- 
ergy (Figure [5]) , which also suggests that the emitting 
area is small. The fitted column density from any model 
in Table Q] agrees reasonably wel l with the T Vh de rived 
from fits to the SNR spectrum. ISun et all (|2004f ) find 
iV H = (1.54 - 1.78) x 10 22 cm" 2 from fitting a vari- 
ety of equilibrium and non-equilibrium ionization mod- 
els to ASCA and C handra spectra of the SNR, while 
IGiacani et all (|2009h find N n = (1.52± 0.02) x 10 22 cm" 2 
from fitting a non-equilibrium ionization model to the 
XMM-Newton observations of 2004. 

We then searched for any indication of variability by 
extracting and fitting the spectrum of each individual 
observation, fixing the spectral parameters (except for 
the total flux) to those of the summed spectrum. In 
order to minimize statistical and systematic errors in this 
comparison, we restricted the fitted energy range to 1 — 
5 keV for the individual spectra. The resulting individual 
fluxes are shown in Figure where they are seen to 
be constant within errors. The mean 1 — 5 keV flux is 
1.9 x 10 -13 erg cm~ 2 s _1 . The typical la uncertainty 
in each observation is w 10%, and the rms dispersion 
among the 23 observations is 10%. 

Unlike the unique case of IE 1207.4-5209, there 
is no evidence for cyclotron absorption lines in the 



spectrum of PSR J1852+0040 or in any other CCO. 
IE 1207.4—5 209 has strong absorption features at 0.7 
and 1.4 keV (jSanwal et alj|2002t iMereghetti et al.ll2002f) 
that can be attributed to the cyclotron fundamental cn- 
ergy E c and its first h armonic in a field of « 8 x 10 10 G 
(Bignami et al.l I2003T ). according to the relation E c = 
1.16(J5/10 n G)/(l + z) keV, where z » 0.3 is the gravi- 
tational redshift. The spectrum of PSR J1852+0040 can 
now be understood in terms of its weaker surface i?-ficld. 
For B s = 3.1 x 10 10 G, the cyclotron fundamental and 
first harmonic are at 0.27 keV and 0.54 keV, below the 
region accessible to X-ray spectroscopy (Figure [3]) , espe- 
cially because of the larger intervening column density 
to PSR J1852+0040. The same explanation applied to 
other well-observed CCOs suggests that they also have 
weaker surface B-fields than IE 1207.4-5209. 

3. DISCUSSION 
3.1. CCOs as Anti-Magnetars 

The steady spin-down of PSR J1852+0040 is consis- 
tent with magnetic braking of an isolated neutron star 
with a weak magnetic field of B s = 3.1 x 10 10 G. In 
this case, the derived characteristic age of 192 Myr is 
not meaningful because the pulsar was born spinning 
at its current period, P = 0.105 s, as was the case 
for the other CCO pulsars, IE 1207.4-5209 and PSR 
J0821-4300, which have P = 0.424 s and P = 0.112 s, 
respecti vely. A recent population analys is of radio pul- 
sars by iFaucher-Giguere fe Kaspl (|2006f ) favors such a 
wide distribution of birth periods (Gaussian centered on 
~ 300 ms, a ~ 150 ms). 

A _B-field of this magnitude could be just the fossil field 
left by flux conservation during the contraction of the 
core of the proge nitor sta r , as or iginally hypothesized for 
neutron stars by IWoltierl (|1964f ). Under this hypothesis, 
CCOs could simply be those neutron stars in which no 
additional mechanism of i?-field amplification has been 
effective. As stronger fields can be generated by a turbu- 
lent dynamo whose stren gth depends on the rotation r ate 
of the proto- neutron star (jThompson fc Duncan|[l993T) . it 
is natural that pulsars born spinning slowly would have 
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Fig. 6. — Fluxes in the 1 — 5 keV band for the 23 individual ob- 
servations of PSR J1852+0040, fitted to the two blackbody model 
given in Table [4] Errors bars are 1 a. The weighted mean flux is 
indicated by the dot-dashed line. 



Spin-down Measurement of PSR J1852+0040 



7 




-22 



A 



0.001 



0.01 



10 



0.1 1 
Period (sec) 

Fig. 7. — P — P diagram of isolated pulsars (dots), binary radio 
pulsars (circled dots), and magnetars (crosses). The location of 
PSR J1852+0040 is marked by the star. The radio pul s ar dea th 
line B/P 2 = 1.7 X 10 11 G s -2 of IBhattacharva et al.l QMS) is 
indicated. The spin-up limit for recycled pulsars corresponds to 
P(ms) = 1.9 (B/10 9 G) 4 / 3 l lvan den HeuvelllT987T) . 



the weaker £>-fields. The model of lBonanno et all (2006) 
supports this, in particular finding that the -B-field of a 
slowly rotating neutron star should be confined to small- 
scale regions, while a global dipole that would be respon- 
sible for spin-down is absent. 

PSR J1852+0040 falls in a region of (P, P) parameter 
space (Figure [7l) that is devoid of ordin ary (non-recycled) 
radio pulsars ( [Manchester et all 120051 ). It overlaps with 
recycled pulsars in this area. The spin parameters of PSR 
J 1852+0040 are not beyond the empirical or theoreti- 
cal radio pulsar death lines of iFaucher-Giguere fc KaspH 
(|2006D or iChen fc Rudermanl (1993D, respectively. One 
possible explanation for the absence of radio emission 
from CCOs is low-level accretion of SN debris for thou- 
sands or even millions of years. However, the sample of 
CCOs is not yet large enough to know if they are intrin- 
sically radio quiet. 

Binary pulsars with similar spin parameters as PSR 
J 1852+0040 are thought to have bee n partly spun up 
by acc retion. It was also suggested bv lDeshpande et all 
(1995) that most single pulsars in this region are recy- 
cled, the birth rate of pulsars with B s < 10 115 G being 
too low in their estimation, perhaps one in ~ 5000 yr. 
To the contrary, IHartman et al.l (|1997f ) concluded that it 
is not possible to distinguish ordinary pulsars of B s = 
l.Qio.5-11.5 from recycled ones, and that there is no rea- 
son to suppose that any such single pulsar is recycled. 
The measurements of spin parameters of CCOs certainly 
increases the inferred birth rate of neutron stars in this 
weak B-field regime, a renewed warning not to assume 
that isolated radio pulsars with similar spin parameters 
are recycled. It should be expected that many of the w 52 
radio pulsars with B s < 10 11 G and 0.1 < P < 0.7 s 
in Figure [7| have characteristic ages that are meaning- 
less, being hundreds of millions of years. They have not 



moved in the (P, P) diagram, and may be former CCOs 
whose SNRs dissipated only 10 5 — 10 6 yr ago. Once the 
SNR has disappeared and the neutron star has cooled, 
it is difficult to recognize and classify such an "orphan 
CCO" using X-rays unless it was first selected as a radio 
pulsar. Allowing for this observational handicap, CCOs 
having weak B-fields may represent a channel of neutron 
star birth that is as common as any other. 

This outcome is counter in its details to the long de- 
bated "injection" hypoth esis. Using the "pulsar cur - 
rent" analysis proposed by iPhinnev fc Blandfordl (1 1981), 
iVivekanand fc Naravanl (|1981[) and lNaravanl (|1987l ) sug- 
gested that a large fraction of pulsars are injected with 
0.5 < P < 1.0 s and P > 1 x 10~ 14 (i.e., with large 
magnetic field), or that neutron stars do not turn on 
as radio pulsars until their P decreases below a critical 
value of ~ 3 x 10~ 13 . This would require that they also 
cool rapidly in order not to be detected as th ermal X-ray 
sources in SNRs. Later i nvestigators (e.g., iLvne et ail 
119851: lLorimer et al.l I1993D argued that this interpreta- 
tion suffers from substantial statistical uncertainties and 
selection effects in pulsar surveys, and that injection of 
a distinct population of pulsars is not necessary. While 
previous authors dismissed the possibilit y of injection of 
radio pulsars with weak B- fields (e.g., [Srinivasa n et all 
1984), this appears to be exactly the origin of radio- 
quiet CCOs. The relat ed problem of t he missing pulsars 
in empty-shell SNRs (jGotthelfl [l998h has largely been 
solved in recent years through a v ariety of channels, as 
summarized by IKaspi fc He lfand (2002), including the 
detection of radio-faint pulsars, pulsar wind nebulae in 
X-rays, and now, pulsations from CCOs. 

The spin-down power E = 3.0 x 10 32 erg s _1 of PSR 
J1852+0040 is an order of magnitude smaller than its 
observed thermal X-ray luminosity, L x « 5.3 x 10 33 dj 1 
erg s _1 , which argues that the latter is mostly residual 
cooling. Remaining questions presented by the X-ray 
properties of PSR J1852+0040 and other CCOs are fo- 
cussed on the details of their X-ray spectra and pulse pro- 
files, which require small heated areas. In the case of PSR 
J1852+0040, the spectrum is more complex than a sin- 
gle blackbody. In the two blackbody model, for example, 
the hotter component, of temperature kT2 = 0.52 keV, 
has an area of only ~ 2.5 dj 1 km 2 . The canonical open- 
field-line polar cap has area A pc — 2tt 2 R 3 /Pc w 1.1 km 2 , 
which is comparable to the area of the hot spot, but po- 
lar cap heating by any magnetospheric accelerator must 
be negligible in this case, being much smaller than the 
spin-down power. Accretion may cover a wider area 
and generate more luminosity, but that process is also 
problematic for PSR J1852+0040, as discussed below. 
In the remainder of this section we discuss alternative 
hypotheses for the weak surface S-fields and thermal 
hot spots of CCOs, involving internal magnetic fields, 
anisotropic conduction, and accretion of fall-back mate- 
rial. Although these theories have some attractive prop- 
erties, they are speculative, and none yet offers a self- 
consistent explanation for all of the observed properties 
of CCOs. 

3.2. Localized Heating by Magnetic Field Decay? 

While the X-ray luminosity of PSR J1852+0040, 5 x 
10 33 erg s _1 , is consistent with minimum neutron star 



Halpern & Gotthelf 



cooling scenarios (iPage et all 12004 [2 0071 for the age of 
Kes 79 (5.4-7.5 kvr. iSun et all 12004 . the high temper- 
atures and small surface areas fitted to the X-ray spec- 
trum are difficult to understand without invoking either 
localized heating on the surface or strongly anisotropic 
conduction. We recall that the X-ray spectra and lu- 
minosities of CCOs are not very different from those of 
quiescent magnetars, another recently recognized class, 
being attributable largely to one or two surface ther mal 
emission components (e.g. JHalpern fc GottheUl2005h . It 
is difficult to distinguish CCOs from quiescent magne- 
tars without timing data or evidence about variability 
(jHalpern fc GotthelJ[2009l ). Thus, it is tempting to hy- 
pothesize that the same magnetic field decay that is 
thought to be responsible for localized crustal heating 
in a magnetar can be operating in CCOs. However, in 
the case of CCOs, such -B-fields would need to have the 
same 10 14 — 10 15 G strengths as in magnetars to account 
for continuous X-ray luminosities of 10 33 — 10 34 erg s" 1 , 
while having insignificant dipole moments that do not 
contribute to spin-down. One such configuration would 
be a small "sunspot" dipole on the rotational equator, 
to account for the large pulse modulation, while the rest 
of the star contributes dipole B s < 3 x f0 10 G in or- 
der not to exceed the observed spin-down rate. While 
we cannot rule out this possibility, it would be remark- 
able if the magnetic field of a neutron star could be cre- 
ated or evolve into a configuration with such extreme 
contrast. Also, the conspicuous lack of X-ray variability 
from CCOs argues against invoking the same magnetic 
field strength and heating mechanism that is held respon- 
sible for magnetars, with their ubiquitous variability. 

3.3. Anisotropic Conduction? 

We next turn to the possibility that magnetic field 
confined entirely beneath the surface can be responsi- 
ble for the nonuniformity of surface temperature. While 
a strong B-field is a favored ingredient of models for 
anisotropic conduction, PSR J1852+0040 with its ex- 
ceptionally weak dipole field would seem the least likely 
candidate for such an explanation. Nevertheless, it is 
possible that a strong toroidal field can exist under the 
surface, while only a weak external poloidal field con- 
tributes to spin-down. A toroidal field is expected to be 
the initial configuration generated by differential rotation 
in the proto-neutron star dynamo ([Thompson fc Duncan! 
f993). The effects on the surface temperature distri- 
bution of an internal tor oi dal field were calculat ed by 
Perez- Azorin et al.l (|2006ft . iGeppert et all ((2006), and 
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field is to insulate the magnetic equator from heat con- 
duction, resulting in warm spots at the poles. It was 
even shown that the warm regions can be of differ- 
ent sizes due to the antisymmetry of the poloidal com- 
ponent of the field, which is reminiscent of the asym- 
metric opposing thermal r egions in the Puppis A CCO 
(IGotthelf fc Halpernl[2009l) . 

In order to explain the large observed pulsed fraction 
of PSR J1852+0040, the axis of toroidal magnetic field 
must make a large angle with respect to the rotation axis 
of the neutron st ar. An orthogonal co nfig uration was 
in fact adopted by IGeppert et al.l (|2006l ) and I Page et al.l 
(2007). Even if the toroidal field component is initially 
parallel to the rotation axis, a toroidal field will deform 



the neutron star into a prolate shape, which tends toward 
orthogonality to the rotation axi s as a minimum energy 
configuration (Braithwaitc 2009). The main problems 
with this mechanism are 1) anisotropic conduction is un- 
likely to produce a small hot region on PSR J1852+0040 
that covers only ~ 1% of the surface area of the neu- 
tron star, and 2) it doesn't explain why only one of the 
orthogonal poles is evidently hot. 

Furthermore, to have a significant effect on the heat 
transport, the crustal toroidal field strength required 
in these models is ~ I0 15 G, many orders of mag- 
nitude greater than the poloidal field if the latter 
is measured by the spin-down. Purely toroidal or 
poloidal fields are thought to be unstable (|Tavlerlll973l ; 
iFlowers fc Rudermanlll977T ), although the toroidal field 
may be stabilized by a polo idal field that is se veral or- 
ders of magnitude smaller (Braithwaitc 200i|), so this 
may be a viable configuration for a CCO. On the other 
hand, twisting and breaking of the crust by such a large 
toroidal field is the basis of the magnetar model for 
soft gamma-ray repeate rs and anomalous X-ray pulsars 
([Thompson et all [2002). which also have large external 
dipole fields as measured by their rapid spin-down. In 
this picture, a CCO is a ma gnetar-in-waiting, a scenario 
that iPavlov fc Lunal (|2009f ) found somewhat contrived, 
as do we. If any pulsar could be an incipient magnetar, 
it isn't explained why CCOs have especially weak surface 
fields compared to ordinary pulsars. 

3.4. Submergence of Magnetic Field by Hypercritical 
Accretion? 

For as long as the SN explosion mechanism has been 
studied, it has been noted that a newly born neu- 
tron star may accrete large amounts of fall-back ma- 
terial in the hours to months after the SN. This could 
be submerge the in itial magnetic field into the core 
(|Geppert et al.|[l999h . and the surface field could be es- 
sentially zero if the accreted matter is not highly mag- 
netized in a well-ordered fashion. After the accretion 
stops ; the submerged field w ill diffuse back to the sur- 
face (JMuslimov & Page 19 951), but this could t ake hun- 
dreds to millions of years ([Geppert et al.|[l999f ). In this 
picture, the CCOs could be those neutron stars that suf- 
fered the most fall-back accretion, while the normal pul- 
sars with ages of a few thousan d years must have a c creted 
<< 0.01 Mq. In the model of iMuslimov fc Page! (|1995f ) 
with accretion of ~ 10~ 5 M Q , the regrowth of the sur- 
face field is largely complete after ~ 10 3 yr, so this would 
probably not explain the weak field of PSR J1852+0040. 
But if > 0.01 Mq is accre ted, then the dif fusion time 
could be millions of years. iChevalier I (|1989f ) calculated 
that the neutron star in SN 1987A could have accreted 
^0.1 Mq of fallback material in the hours after the SN 
explosion, aided by a reverse shock from the helium layer 
of the progenitor. If so, it may never emerge as a ra- 
dio pulsar. A more normal type II SN should accrete 
10~ 3 Mq or less, which may delay the emergence of the 
magnetic field for tens of thousands of years, a timescale 
applicable to explaining the properties of CCOs. 

Finally, it should be considered that a thermoelectric 
instability mechanism drive n by the strong tempe rature 
gradient in the outer crust fBlandfo rd et al.| [T983) could 
regenerate the magnetic field submerged by accretion. 
However, the growth of the field via this mechanism could 
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take ~ 10 5 yr. CCOs could then be those neutron stars 
in which the thermoelectric instability is the dominant 
mechanism of generating magnetic field, perhaps because 
their initial rotation rates are slow. 

The models discussed here have the effect of delaying 
the emergence of a normal magnetic field for times rang- 
ing from a thousand years to essentially forever. They 
arc difficult to test using CCOs, for which it is impossible 
to measure the braking index, which could indicate field 
amplification. Furthermore, they don't help to explain 
the small, hot surface areas seen in the X-ray spectra 
of CCOs. In order to achieve such a configuration, the 
magnetic field would have to be submerged everywhere 
except for one or two small regions, perhaps the mag- 
netic poles of the neutron star, which could remain hot 
via conduction from the interior. 

3.5. Continuing Accretion from a Fall-Back Disk? 

We also consider ongoing accretion from a fall- 
back disk of supernova debris as a possible source of 
anisotropic surface heating. While variability is an in- 
dicator of accretion, we do not have evidence of any vari- 
ations of PSR J1852+0040 at the 10% level among 23 
observations spanning 4.8 yr, which would tend to elimi- 
nate accretion as a significant contributor to its X-ray lu- 
minosity. However, although accretion is widely consid- 
ered to be an inherently unstable process, it is not clear 
that variability of X-ray binaries or AGNs can be ex- 
trapolated to accretion from a fossil disk at rates of only 
~ 10~ 5 ^Edd- Therefore, we also explore the possibility 
of accretion, constrained only by the steady spin-down 
rate of PSR J1852+0040, using the theory of propeller 
and accretion-disk torques. 

The spin parameters of PSR J1852+0040 fall in a 
regime in which both dipole braking and accretion 
torques are conceivably significant. In Paper 2, we dis- 
cussed the implications for P of PSR J1852+0040 accret- 
ing and spinning down in the propeller regime, deriving 
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for the propeller effect. In this model, M is the rate of 
mass expelled, which must be > m, which is the accre- 
tion rate onto the neutron star. We can therefore assume 
Lx = rjrnc 2 to calculate a lower limit on M in the pro- 
peller accretion scenario. Assuming efficiency 77 ~ 0.1, 
M > 5 x 10 13 g s _1 is required. But if we were to adopt 

= B S R 3 , where B s = 3.1 x 10 10 G from assuming 
dipole spin-down, then the observed P — 8.68 x 10~ 18 
allows a negligible M < 3 x 10 8 g s _1 , which contradicts 
the accretion assumption. 

If instead we use the measured P as an upper limit on 
the propeller spindown rate, we can derive an upper limit 
on the required \x in the accretion scenario. Under these 
assumptions, ^< 3.4x 10 26 G cm 3 , or B s < 3.4 x 10 9 G. 
However, this would marginally violate the conditions of 
the propeller model, because for such a small magnetic 
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field, the magnetospheric radius is 

r m = 2.3xl0 7 4 7 Mr 3 2/7 
which is smaller than the corotation radius 

r co = (GM) 1/3 (^j ' = 3.7 x 10 7 cm, 

and the pulsar would be rotating near its equilibrium pe- 
riod. But there has not been enough time in the life of 
Kes 79 for a weakly accreting pulsar to come into equilib- 
rium. To find PSR J1852+0040 in an equilibrium state 
would be a remarkable coincidence because its spin-down 
timescale, P/P ~ 4 x 10 8 yr, is orders of magnitude 
greater than its age. Therefore, it is more natural to 
conclude that the small P of PSR J1852+0040 is due 
purely to dipole spin-down, and that it is not accreting. 
Absence of accretion would also argues against nonuni- 
form surface composition as a cause of the temperature 
variations. Finally then, none of the possible mechanisms 
discussed here for hot spots on CCOs is entirely natural 
and self-consistent. 

4. CONCLUSIONS AND FURTHER WORK 

Measurement of the spin-down rate of the 105 ms PSR 
J 1852+0040 in Kes 79 was achieved using X-ray obser- 
vations coordinated between XMM-Newton and Chan- 
dra. The resulting phase-connected ephemeris span- 
ning 4.8 yr requires only the first derivatives of the 
spin frequency, with no measurable timing noise super- 
posed. The straightforward interpretation of this result 
is dipole spin-down due to a weak surface magnetic field 
of B s = 3.1 x 10 10 G, the smallest measured for any young 
neutron star. While this is the first measurement of the 
spin-down rate of a CCO, upper limits on the period 
derivatives of IE 1207.4-5209 and PSR J0821-4300, as 
well as spectral features in those two, also indicate weak 
S-fields. The properties of the CCO class were loosely 
defined in the past, but it is now evident that a weak 
S-field is the physical reason that a large fraction of the 
neutron stars were so named. The body of evidence sup- 
ports the "anti-magnctar" model for the origin of such 
CCOs: neutron stars born spinning "slowly" may as a 
result be endowed with only weak magnetic fields. PSR 
J 1852+0040 falls in a region of (B, P) space that over- 
laps with moderately recycled pulsars. Single radio pul- 
sars with similar spin parameters may therefore be for- 
mer CCOs rather than recycled, and they may be much 
younger than their characteristic ages. X-ray observa- 
tions of pulsars in this region may find evidence of their 
relative youth via surface thermal emission. 

Ongoing timing studies of the CCOs IE 1207.4-5209 
and PSR J0821— 4300 will able to measure their spin- 
down rates and can refine the correspondence be- 
tween spectral features interpreted as electron cyclotron 
lines and surface 5-ficld. Those CCOs such as PSR 
J1852+0040, with no apparent absorption lines, may 
simply have weaker £?-fields than the others by a factor 
2 — 3, so that the cyclotron resonance is below the soft 
X-ray band. As for CCOs that have not yet been seen 
to pulse, there is no reason to suppose that they are of 
a different physical class. Rather, their surface tempera- 
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ture distributions and/or viewing geometries may deliver 
only weak modulation that happens to fall below the sen- 
sitivity of existing data, while deeper observations may 
succeed in discovering their spins. It would be of great 
interest to discover pulsations from the CCO in Cas A, 
the youngest known neutron star, which in all respects is 
the prototype of the CCO class. At an age of 330 yr, it 
is an order of magnitude younger than the others, whil e 
having similar spectral properties (|Pavlov fc Lunafl 2009) . 
In our model, its spin parameters were fixed at birth, and 
should conform closely to those of the older CCOs. A re- 
cent Chandra observation of the Cas A CCO did not de- 
tect pulsations (see Appendix), but the High Resolution 
Camera (HRC) that was employed lacks energy resolu- 
tion, a capability that may be crucial in the same way 
that it was fo r the XMM-Newton discov ery of the pulsar 
in Puppis A (jGotthelf fe Halperril2009l) . 

A remaining theoretical puzzle about CCOs is the ori- 
gin of their surface temperature anisotropies, in partic- 
ular, the one or two warm/hot regions that are smaller 
than the full neutron star surface. The high signal-to- 
noise spectrum accumulated from PSR J1852+0040 re- 
veals similar temperature structure as the other CCOs. 
The measured spin-down power of PSR J1852+0040 is 
too small by an order of magnitude to contribute to these 
excess emissions. We considered whether accretion of 
fall-back material can be responsible for this effect, but 
find it unlikely in the case of PSR J1852+0040 because 
of its slow, steady spin-down. If other CCOs are found to 
have similar spin-down properties, the same arguments 
would apply to them. Therefore, it is important to mea- 
sure the spin-down rate of PSR J0821-4300, the CCO 
in Puppis A, and investigate furthe r its apparent phase- 
depen dent emission line at 0.8 keV (jGotthelf fc Halpern! 
l20"09f) . which may be indicative of accretion because it is 
in emission. 

The small regions of high surface temperature are prop- 



erties that CCOs share with magnetars, although most 
magnetars are hotter and more luminous. The total X- 
ray luminosities of CCOs are consistent with slow cooling 
scenarios, and there would be no need to hypothesize a 
mechanism of magnetic field amplification beyond simple 
flux freezing if it were not for their strongly anisotropic 
surface temperature distributions. Paradoxically, the ex- 
planations that have been considered for hot spots invoke 
strong magnetic fields just below the surface, which is 
exactly the basis of the magnetar model. In this sense, 
CCOs are potentially magnetars waiting to emerge, or 
maybe ordinary neutron stars in which the original field 
was submerged by especially massive fall-back of super- 
nova debris. This is an unsatisfying picture, if only be- 
cause the contrast between field strength in different re- 
gions must be a factor of 10 4 or larger to affect the sur- 
face temperature distribution while maintaining the slow 
spin-down and allowing soft X-ray cyclotron lines. Re- 
solving this essential mystery of the CCOs will undoubt- 
edly generate new insight into the physics of neutron 
stars. 
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APPENDIX 

SEARCH FOR X-RAY PULSATIONS FROM THE CAS A CCO 

In this Appendix we report on analysis of a recent Chandra observation designed to be the most sensitive yet to 
search for pulsations from CXO U J232327. 9+584842, the CCO in Cas A. But first we summarize the results of previous 
attempts. iMurrav et al.l (|2002f ) searched a 50 ks Chand ra HRC observ a tion and identified sev eral candidate periods, 
but none was confirmed in a follow-up 50 ks observation ((Ransom 2002T). iPavlov fc Lu na (2009) found from a Chandra 
ACIS observation taken in subarray mode / p (3ct) < 16% for P > 0.7 s. That mode is limited by its sampling time of 
0.34 s. ACIS in full-frame mode has a sampling time of 3.2 s, so is only sensitive to periods > 6.4 s. For the 1 Ms 
of data taken in this mode ([Huang eTaDHQOl we searched the range 6.4 < P < 100 s and / < 1 x 10 11 Hz s 1 , 
corresponding to t c > 330 yr. The resulting upper limit is f p < 4% at 99% confidence. Existing XMM-Newton data 
on CXOU J232327. 9+584842 are also limited in sensitivity and time resolution, and do not cover the full range of 
period s spanned by known C CO pulsars. Using an effective exposure of 9.2 ks in the EPIC pn taken in full-frame 
mode, iMereghetti et al.l (|2002D reported an upper limit of f p < 13 % for P > 0.3 s, and f v < 7% for P > 3 s. However, 
these limits may have been underestimated by a factor of 2 — 3 ( see IPavlov fc Lunall2009f ). Our own analysis sets limits 
of f p < 32% for 0.146 ms < P < 3 s and f p < 17% for P > 3 s from this observation. Improvements in sensitivity and 
time resolution could be obtained with a longer XMM-Newton observation using the EPIC pn in SW mode, as was 
employed for the other CCO pulsars. 

Recently, 487 ks of Chandra HRC timing data were obtained over a period of 11 days in 2009 March, of which 433 ks 
are in the public archive. We searched the public data for pulsations. A total of 11,486 counts were extracted from 
a circular aperture of radius l."2, of which 6% are estimated to be background. We sampled (/, /) parameter space 
using the Zf test, which is an optimal one for a source in which the light curve, presumed to arise from surface thermal 
emission, is expected to be weakly modulated and quasi-sinusoidal. We limited the search range to / < 200 Hz for 
/ < 5 x 10~ 13 Hz s - \ / < 100 Hz for / < 3 x 10~ 12 Hz s - \ and / < 10 Hz for / < 1 x 10~ 10 Hz s" 1 , oversampling 
by a factor of « 3 in each parameter. This range covers the parameters of all the known CCO pulsars and magnetars. 
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The largest values of Z\ found were ~ 40. The theoretical distribution of noise power is exponential with a mean of 
2, so the single trial probability that Z\ > 40 by chance is 2 x 10~ 9 . Such a value is expected to arise randomly in 
the ~2x 10 9 independent trials that were performed, so it is not significant. The expectation value of the intrinsic 
pulsed fraction is f p = (1 + N b /N s )^/2ZjjN, where N b /N s = 0.06/0.94 is the ratio of background to source counts, 
and N = N s + Ni, = 11,486 is the total number of counts. In this case, a sinusoidal signal of f p — 9.0% would have 
a 50% probability of giving Z\ > 40. However, taking into account the noise, which may either increase or decrease 
the total power, the 99% upper limit on the Z\ of a true signal in this case is 75 <|Grothlll975t IVaughan et al1fl99l . 
corresponding to f p < 12.2%. On the other hand, these peak power values are somewhat overestimated because we 
oversampled Fourier space. Finally, then, we adopt an upper limit of f p < 12% for P > 10 ms in Table [1] which is the 
lowest limit yet obtained for Cas A over the range of periods that is spanned by the known CCO pulsars. 
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